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Abstract 

The creation and annihilation operators of a two-term diatomic molecular potential are studied 
and it is observed that they satisfy the commutation relations of a SU(1,1) algebra. To study the 
Lie algebraic realization of the present potential, the normalized eigenfunctions and eigenvalues 
are computed by using the Nikiforov-Uvarov method. 
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I. 



INTRODUCTION 



The two-term diatomic molecular potential under consideration [1] 



V{r) = -Vo 




(1) 



1 - qe-^' 



(1 - qe-^^) 



2 ' 



has been firstly proposed by Sun to study of some diatomic molecules. Jia and co-workers 
have studied the bound states of this potential for the Schrodinger (SE) equation in the 
view of supersymmetric approach [2] and in Ref . [3] , the energy spectra of the potential is 
obtained in terms of Green's function for the same domain. Furthermore, Arda et. al. have 
studied the bound states and the corresponding normalized wave functions of this potential 
for the SE equation with any ^-values for the cases where q > 1 and q = 0. In addition, 
they have also given the analytical results for the Manning-Rosen potential, the Hulthen 
potential and the generalized Morse potential as special cases with their numerical results 
[4] . As a continuation, in the present work, we intend to search the creation and annihilation 
operators of the two-term molecular potential and find that the dynamical group satisfied by 
the ladder operators of the present potential. For this aim, we obtain firstly the normalized 
eigenfunctions of the potential by using the Nikiforov-Uvarov (NU) method [5] to present 
the ladder operators. 

The group theoretical realizations within the quantum mechanics have been received great 
attentions in literature. The factorization method is used to get exact solutions of some 
solvable potentials with the ladder operators [6, 7]. The perturbed ladder operator method 
is applied to quantum mechanical systems with the help of unperturbed eigenvalues and 
eigenfunctions [6, 8]. Dong and co-workers have applied a 'different' factorization method 
to various types of molecules [9-16]. 

The organization of the present work is as follows. In Section 2, we search the normalized 
wave functions with their energy spectra of this potential by using NU-method which is 
briefly given in Appendix A. In Section 3, we present the raising and lowering operators of 
the potential and show that the commutators of the ladder operators satisfy the Lie algebra 
of a SU(1,1) group. 
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II. NORMALIZED WAVE FUNCTIONS AND ENERGY SPECTRUM 



The radial Schrodinger equation is written [17] 



2/1 £{£ + !) 



[Ene - V{r)] - 



Rne{r) = 0, (2) 



where i is the angular momentum quantum number, m is the particle mass, V{r) is the 
central potential and Ene is the non-relativistic energy. 

Inserting Eq. (1) into Eq. (2) and using a new variable x = qe"^'^ (x — )■ g for r — >■ 0, 
a; ^ for r — )■ oo) gives the following equation without the centrifugal term 

d Rn{x) 1 — X dRnix) r 2 2 2 2^ id t \ n 

+ + - - ] i?n(x) = . (3) 

with the abbreviations 

% = , (4a) 



a^ = -^ 2K-^ , (4b) 



2 _ 2/x ( ^ ^ Vo 

"-J^(^7-")' 

where Vq, Vi, ^ and q are real parameters defined by Vi = Do{e^ ~ q): Vo — 2Vi, ^ = A*/ro, 
where Dq is the depth of the potential, Tq is the equilibrium of the molecule and q is the 
shape parameter. 

Comparing Eq. (3) with Eq. (Al)we obtain 

f{x) = 1 — X , (5a) 
a{x) = X ~ x'^ , (5b) 
a{x) — —a^x^ — a\x — a\ , (5c) 

with a"{x) = -2. Eq. (A7) gives 



n{x) = =F y + a| - /c^ a;2 + (a| + /c) x + af , (6) 

We obtain the parameter k is determined by setting the discriminant of the square root in 
Eq. (6) 

ki — —al — 2a\ — aiA , (7a) 



k2 = -al-2al + aiA- ^ = J ^ + • C^^) 
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Prom Eqs. (7a) and (7b) we get the foUowings for 7r(x), respectively, 

7r[x) = ai — [ Oi H — ) x , (8aj 

(8b) 

The parameter A required for the method is obtained from Eq. (AS) with the help of Eq. 
8(a) 



7r(a;) = ai — I ai H — l x . 



A = -al - 2al - {ai + -] - A{ai + -] , 



and the other one A„ is written 



An = n (2 + 2ai + A) + n(n - 1) , 



(9) 



(10) 



where used t{x) = 1 + 2ai - (2 + 2ai + A)x and t'{x) = -(2 + 2ai + A). 

Setting A = A„ and using the parameter values given in Eqs. (4a)- (4c) gives the following 
energy eigenvalue equation 



ai + 



- in + - J ^ 03 , 



(11) 



and the analytical expression for the energy spectra 



En — 



2ii 





^ 8m (Va ,Vi\ 
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2n + l + ^l + ^_ 





(12) 



In order to find the wave functions, we first use Eq. (A6) to present the function p{x) 

p(x) ~x2"i(l-x)^, (13) 
and inserting it into Eq. (A5) gives 



(14) 



a;2«i(l - x)^ dx" 

which means that the second part of the wave function in Eq. (A2) is written as the Jacobi 
polynomials [18] 



^n{x)^Pt'^^\l-2x), 
Eq. (A4) gives the first part of wave function as 

^(x) ~x'^^(l-x)(^+^)/2^ 
4 



(15) 



(16) 



Finally, the total wave functions are given as 

Rn{x) = AnX^'il - x)(^+^)/2pi2"i'^)(l - 2x) . (17) 

where An is a constant determined from the normalization condition. Choosing a new 
variable as x — qy, we write the normalization condition 

\Anf q^'^^y^'^^il - qyf^^ [P^^a^.A)^! _ 3^^)] [p(^-^^^){l - 2qy)\ dy = 1 , (18) 

^0 

Using the following representation of the Jacobi polynomials [18] 

,c ^ . 1 ( a + n \ I b + n \ , , 

= { ]{i+zr{i-zr-^ (19) 

^ k=o \ k J \n-kj 
where binomial coefficient ^ ^ = ^K^^-ry. ~ r(r+i)r(n-r+i) ' Hence, Eq. (18) becomes 

\Anf g2ai+n+m-fe-^+l^_^y«+m-fc-£ ^^^^^ ^ ^^^^ f ^2oi+n-fc+m-£^^ _ qyf+A+k+t^y ^ ^ ^ 

Jo 

(20) 

where g{n, k) and 5'(m, t) are two arbitrary functions of the parameters Oi and A and given 

by 

( n + 2ai \ ( n + A\ 
fe=o \ J \n K J 

g{m,i) ^ g{n,k){n ^ m;k ^ i) , (21b) 

By using the following integral representation of hypergeometric type function 2P1 {a, b; c; z) 
[18] 



2Fi(a, 6; c; ^) = t'"^ (1 - t)^"''"^ (1 - tz)-'^ dt , 



(22) 



and by setting the variable z ^ q and taking a = l + A + k + i, b = 2ai + n + 'm — k — £+l, 
c — 1 + b, Eq. (20) gives the normalization constant as 



^ _ ,'g-2»i-"-"^+'=+^-i(-l)fc+^-"-™r(2ai+n + m-A;-£ + 2) 

5(n,A;)5(m,£)r(2ai + n + m- A;-£+l)2Pi(a,6;c;g) ^ 



III. CREATION AND ANNIHILATION OPERATORS 



We want to obtain the following eigenvalue equations for creation, fi+, and annihilation, 
n_, operators 



(24) 



where the ket \n > corresponds to the eigenfunctions Rn{x). Firstly, we study the effect of 
the differential operator ^ on the eigenfunctions 
d 



dx 



\n >-- 



With the help of the following relation [18] 



I. [pM{x)] = l(n + rj, + rj, + l)Pi'2^'''^'-'\x) 



we obtain the first derivative in the last term of Eq. (25) as 



dx 



[P(2"-^)(1 - 2a;)] = -(n + 2a,+A+ 1)pS 



(2ai+l,A+l) 



Xj 



1 



2n + rji + r]2 + 2 1- x 



1 



Using the recursion relations for the Jacobi polynomials [18] 
^ ' m + n + l " ^ ^ ryi+n + 1" 

Eq. (25) becomes 

d 



2n + ?7i + ?72 + 2 1 + X 



n 



(n 



+ V2 + l)Pi'''''^\x) + {n + l)Pi'T\x) 



dx 



\n >= 



l+A 

dl - 1^3 2 '^2 



X 



1-X 



n > + — + 



Kl Ki \ A. 



with 



{n + A){n + 2ai) 



n(n + 2a 



X 1 X J -^ji 1 

n(n + A) 



n-l> 



2n + 2ai + A 2n + 2ai + A 

Rewriting Eq. (29) in the standard form 



2n + 2ai + A 



(25) 



(26) 



(27) 



(28a) 



(28b) 



(28c) 



(29) 



(30) 



l + A 



- K2 - K,3 ] X - ai + K3 



,2n + 2ai + A., , An , 

X )\n>^ n + AH^|n-l>, 

n + zoi An-i 



(31) 
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we obtain the annihilation operator 

n_ 

with the eigenvalues 



d 1 

x{l - x)— + - (2ai + A - 2{k2 + K3) + l)x - ai + K3 
ax 2 



,2n + 2ai + A 
' n + 2ai 



n-^{n + A)-^. 

which satisfy the eigenvalue equation in Eq. (24). 

Following the same procedure and using the following relation [18] 



), (32) 



(33) 



X] 



we immediately write the creation operator as 



X{1 ~ + ^ (2^1 + ^ + 2(^4 - ^5) + 1)2; - Oi - /t4 



,2n + 2ai + A + 2 



where 



_ n + 2ai + A + l)(n + 2ai + 1) 



2n + 2ai + A + 2 
with the eigenvalues 



n + 2ai + A+l 

n + 2ai + A + l)(n + A + 1) 
2n + 2ai + A + 2 ' 



7r+ = (n + 1) 



(34) 



) , (35) 



(36) 



(37) 



Now, let us study the commutation relations of the operators Hip to present the Lie algebra 
related with the two-term molecular potential. For this aim, we define the operator 

1 + A 



Ho = n + 



where n is the number operator satisfying 



2 ' 



n\n >— n\n > , 



So, the operator Ho has the eigenvalues 



ttq = n + 



1 + A 



(38) 



(39) 



(40) 



By using Eqs. (32), (35), (38) and with the help of Eqs. (33), (37) and (40), we obtain the 
commutators as 



n_,n+ 



= 2no ; 



n_,no 



(41) 
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which correspond to the Lie algebra associated with a SU(1,1) group. 

Finally, we could write the Casimir operator of the group by using the above results as 

c = -n_fi+ + Ho (Ao + 1) = -n+n- + no (Ao - 1) , (42) 

with the eigenvalue equation 

C\n >= Co\n > ; (43) 

where 

1 + A 

Co = c(c - 1) ; c = . (44) 



IV. CONCLUSIONS 



We have computed the exact bound state energy eigenvalues and the corresponding nor- 
malized eigenfunctions of the two-term diatomic potential by using the NU method. Then 

we have obtained the step-up and step-down operators of the potential and shown that the 
ladder operators satisfy the Lie algebra of a SU(1,1) group. The Casimir operator of the 
group is also obtained and the eigenvalue equation is given. 
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Appendix A: Nikiforov Uvarov Method 

The Schrodinger equation can be transformed by using appropriate coordinate transfor- 
mation into following form 

where a{y) and d-{y) arc polynomials, at most, second degree, and f{y) is a first-degree 
polynomial. By using the separation of variables, the solution is written as 

*(?/) = 'iPiyMy) , (A2) 



which gives Eq. (Al) as a hypergeometric type equation [5] 

c^V(z/) , r{y)d(p{y) , A 



where ip{y) is defined by using the equation [5] 

Hv) dy a{y) ' 
and the other part of the solution in Eq. (A2) is given by 

where a„ is a normahzation constant and p{y) is the weight function and satisfies the fol- 
lowing equation [5] 

Hy)p{y)]' = r{y)p{y) . (A6) 

where prime denotes the derivative to coordinate. 

The function Tr{y) and the parameter A in the above equation are defined as 



Ay) = I [Ay) - r{y)] ± ^J\ W'{y)-f{y)?-Ay) + ka{y) , (A7) 

A = A; + 7r'(y) . (A8) 

In the NU method, the square root in Eq. (A7) must be the square of a polynomial, so the 
parameter k can be determined. Thus, a new eigenvalue equation becomes 

A = A„ = -nr'{y) -\n{n- \)a"{y) . (A9) 

where prime denotes the derivative and the derivative of the function T(y) = f(y) + 27r(y) 
should be negative. 
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